\(\int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx\) [3]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [A] (verification not implemented)
   Sympy [B] (verification not implemented)
   Maxima [A] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 30, antiderivative size = 77 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=\frac {1}{8} a^2 c x-\frac {a^2 c \cos ^3(e+f x)}{3 f}-\frac {a^2 c \cos (e+f x) \sin (e+f x)}{8 f}+\frac {a^2 c \cos (e+f x) \sin ^3(e+f x)}{4 f} \]

[Out]

1/8*a^2*c*x-1/3*a^2*c*cos(f*x+e)^3/f-1/8*a^2*c*cos(f*x+e)*sin(f*x+e)/f+1/4*a^2*c*cos(f*x+e)*sin(f*x+e)^3/f

Rubi [A] (verified)

Time = 0.08 (sec) , antiderivative size = 77, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.167, Rules used = {3045, 2718, 2715, 8, 2713} \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=-\frac {a^2 c \cos ^3(e+f x)}{3 f}+\frac {a^2 c \sin ^3(e+f x) \cos (e+f x)}{4 f}-\frac {a^2 c \sin (e+f x) \cos (e+f x)}{8 f}+\frac {1}{8} a^2 c x \]

[In]

Int[Sin[e + f*x]*(a + a*Sin[e + f*x])^2*(c - c*Sin[e + f*x]),x]

[Out]

(a^2*c*x)/8 - (a^2*c*Cos[e + f*x]^3)/(3*f) - (a^2*c*Cos[e + f*x]*Sin[e + f*x])/(8*f) + (a^2*c*Cos[e + f*x]*Sin
[e + f*x]^3)/(4*f)

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2713

Int[sin[(c_.) + (d_.)*(x_)]^(n_), x_Symbol] :> Dist[-d^(-1), Subst[Int[Expand[(1 - x^2)^((n - 1)/2), x], x], x
, Cos[c + d*x]], x] /; FreeQ[{c, d}, x] && IGtQ[(n - 1)/2, 0]

Rule 2715

Int[((b_.)*sin[(c_.) + (d_.)*(x_)])^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((b*Sin[c + d*x])^(n - 1)/(d*n))
, x] + Dist[b^2*((n - 1)/n), Int[(b*Sin[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n, 1] && Integ
erQ[2*n]

Rule 2718

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3045

Int[sin[(e_.) + (f_.)*(x_)]^(n_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + (f_.
)*(x_)]), x_Symbol] :> Int[ExpandTrig[sin[e + f*x]^n*(a + b*sin[e + f*x])^m*(A + B*sin[e + f*x]), x], x] /; Fr
eeQ[{a, b, e, f, A, B}, x] && EqQ[A*b + a*B, 0] && EqQ[a^2 - b^2, 0] && IntegerQ[m] && IntegerQ[n]

Rubi steps \begin{align*} \text {integral}& = \int \left (a^2 c \sin (e+f x)+a^2 c \sin ^2(e+f x)-a^2 c \sin ^3(e+f x)-a^2 c \sin ^4(e+f x)\right ) \, dx \\ & = \left (a^2 c\right ) \int \sin (e+f x) \, dx+\left (a^2 c\right ) \int \sin ^2(e+f x) \, dx-\left (a^2 c\right ) \int \sin ^3(e+f x) \, dx-\left (a^2 c\right ) \int \sin ^4(e+f x) \, dx \\ & = -\frac {a^2 c \cos (e+f x)}{f}-\frac {a^2 c \cos (e+f x) \sin (e+f x)}{2 f}+\frac {a^2 c \cos (e+f x) \sin ^3(e+f x)}{4 f}+\frac {1}{2} \left (a^2 c\right ) \int 1 \, dx-\frac {1}{4} \left (3 a^2 c\right ) \int \sin ^2(e+f x) \, dx+\frac {\left (a^2 c\right ) \text {Subst}\left (\int \left (1-x^2\right ) \, dx,x,\cos (e+f x)\right )}{f} \\ & = \frac {1}{2} a^2 c x-\frac {a^2 c \cos ^3(e+f x)}{3 f}-\frac {a^2 c \cos (e+f x) \sin (e+f x)}{8 f}+\frac {a^2 c \cos (e+f x) \sin ^3(e+f x)}{4 f}-\frac {1}{8} \left (3 a^2 c\right ) \int 1 \, dx \\ & = \frac {1}{8} a^2 c x-\frac {a^2 c \cos ^3(e+f x)}{3 f}-\frac {a^2 c \cos (e+f x) \sin (e+f x)}{8 f}+\frac {a^2 c \cos (e+f x) \sin ^3(e+f x)}{4 f} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.11 (sec) , antiderivative size = 47, normalized size of antiderivative = 0.61 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=\frac {a^2 c (12 e+12 f x-24 \cos (e+f x)-8 \cos (3 (e+f x))-3 \sin (4 (e+f x)))}{96 f} \]

[In]

Integrate[Sin[e + f*x]*(a + a*Sin[e + f*x])^2*(c - c*Sin[e + f*x]),x]

[Out]

(a^2*c*(12*e + 12*f*x - 24*Cos[e + f*x] - 8*Cos[3*(e + f*x)] - 3*Sin[4*(e + f*x)]))/(96*f)

Maple [A] (verified)

Time = 0.79 (sec) , antiderivative size = 46, normalized size of antiderivative = 0.60

method result size
parallelrisch \(-\frac {a^{2} c \left (-12 f x +24 \cos \left (f x +e \right )+3 \sin \left (4 f x +4 e \right )+8 \cos \left (3 f x +3 e \right )+32\right )}{96 f}\) \(46\)
risch \(\frac {a^{2} c x}{8}-\frac {a^{2} c \cos \left (f x +e \right )}{4 f}-\frac {a^{2} c \sin \left (4 f x +4 e \right )}{32 f}-\frac {a^{2} c \cos \left (3 f x +3 e \right )}{12 f}\) \(60\)
derivativedivides \(\frac {-a^{2} c \left (-\frac {\left (\sin ^{3}\left (f x +e \right )+\frac {3 \sin \left (f x +e \right )}{2}\right ) \cos \left (f x +e \right )}{4}+\frac {3 f x}{8}+\frac {3 e}{8}\right )+\frac {a^{2} c \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )}{3}+a^{2} c \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-a^{2} c \cos \left (f x +e \right )}{f}\) \(106\)
default \(\frac {-a^{2} c \left (-\frac {\left (\sin ^{3}\left (f x +e \right )+\frac {3 \sin \left (f x +e \right )}{2}\right ) \cos \left (f x +e \right )}{4}+\frac {3 f x}{8}+\frac {3 e}{8}\right )+\frac {a^{2} c \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )}{3}+a^{2} c \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )-a^{2} c \cos \left (f x +e \right )}{f}\) \(106\)
parts \(-\frac {a^{2} c \cos \left (f x +e \right )}{f}+\frac {a^{2} c \left (-\frac {\cos \left (f x +e \right ) \sin \left (f x +e \right )}{2}+\frac {f x}{2}+\frac {e}{2}\right )}{f}+\frac {a^{2} c \left (2+\sin ^{2}\left (f x +e \right )\right ) \cos \left (f x +e \right )}{3 f}-\frac {a^{2} c \left (-\frac {\left (\sin ^{3}\left (f x +e \right )+\frac {3 \sin \left (f x +e \right )}{2}\right ) \cos \left (f x +e \right )}{4}+\frac {3 f x}{8}+\frac {3 e}{8}\right )}{f}\) \(114\)
norman \(\frac {-\frac {2 a^{2} c}{3 f}+\frac {a^{2} c x}{8}-\frac {2 a^{2} c \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 f}-\frac {2 a^{2} c \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{f}-\frac {2 a^{2} c \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{f}-\frac {a^{2} c \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{4 f}+\frac {7 a^{2} c \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}-\frac {7 a^{2} c \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}+\frac {a^{2} c \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4 f}+\frac {a^{2} c x \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}+\frac {3 a^{2} c x \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{4}+\frac {a^{2} c x \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{2}+\frac {a^{2} c x \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{8}}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{4}}\) \(244\)

[In]

int(sin(f*x+e)*(a+a*sin(f*x+e))^2*(c-c*sin(f*x+e)),x,method=_RETURNVERBOSE)

[Out]

-1/96*a^2*c*(-12*f*x+24*cos(f*x+e)+3*sin(4*f*x+4*e)+8*cos(3*f*x+3*e)+32)/f

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 63, normalized size of antiderivative = 0.82 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=-\frac {8 \, a^{2} c \cos \left (f x + e\right )^{3} - 3 \, a^{2} c f x + 3 \, {\left (2 \, a^{2} c \cos \left (f x + e\right )^{3} - a^{2} c \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{24 \, f} \]

[In]

integrate(sin(f*x+e)*(a+a*sin(f*x+e))^2*(c-c*sin(f*x+e)),x, algorithm="fricas")

[Out]

-1/24*(8*a^2*c*cos(f*x + e)^3 - 3*a^2*c*f*x + 3*(2*a^2*c*cos(f*x + e)^3 - a^2*c*cos(f*x + e))*sin(f*x + e))/f

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 245 vs. \(2 (70) = 140\).

Time = 0.18 (sec) , antiderivative size = 245, normalized size of antiderivative = 3.18 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=\begin {cases} - \frac {3 a^{2} c x \sin ^{4}{\left (e + f x \right )}}{8} - \frac {3 a^{2} c x \sin ^{2}{\left (e + f x \right )} \cos ^{2}{\left (e + f x \right )}}{4} + \frac {a^{2} c x \sin ^{2}{\left (e + f x \right )}}{2} - \frac {3 a^{2} c x \cos ^{4}{\left (e + f x \right )}}{8} + \frac {a^{2} c x \cos ^{2}{\left (e + f x \right )}}{2} + \frac {5 a^{2} c \sin ^{3}{\left (e + f x \right )} \cos {\left (e + f x \right )}}{8 f} + \frac {a^{2} c \sin ^{2}{\left (e + f x \right )} \cos {\left (e + f x \right )}}{f} + \frac {3 a^{2} c \sin {\left (e + f x \right )} \cos ^{3}{\left (e + f x \right )}}{8 f} - \frac {a^{2} c \sin {\left (e + f x \right )} \cos {\left (e + f x \right )}}{2 f} + \frac {2 a^{2} c \cos ^{3}{\left (e + f x \right )}}{3 f} - \frac {a^{2} c \cos {\left (e + f x \right )}}{f} & \text {for}\: f \neq 0 \\x \left (a \sin {\left (e \right )} + a\right )^{2} \left (- c \sin {\left (e \right )} + c\right ) \sin {\left (e \right )} & \text {otherwise} \end {cases} \]

[In]

integrate(sin(f*x+e)*(a+a*sin(f*x+e))**2*(c-c*sin(f*x+e)),x)

[Out]

Piecewise((-3*a**2*c*x*sin(e + f*x)**4/8 - 3*a**2*c*x*sin(e + f*x)**2*cos(e + f*x)**2/4 + a**2*c*x*sin(e + f*x
)**2/2 - 3*a**2*c*x*cos(e + f*x)**4/8 + a**2*c*x*cos(e + f*x)**2/2 + 5*a**2*c*sin(e + f*x)**3*cos(e + f*x)/(8*
f) + a**2*c*sin(e + f*x)**2*cos(e + f*x)/f + 3*a**2*c*sin(e + f*x)*cos(e + f*x)**3/(8*f) - a**2*c*sin(e + f*x)
*cos(e + f*x)/(2*f) + 2*a**2*c*cos(e + f*x)**3/(3*f) - a**2*c*cos(e + f*x)/f, Ne(f, 0)), (x*(a*sin(e) + a)**2*
(-c*sin(e) + c)*sin(e), True))

Maxima [A] (verification not implemented)

none

Time = 0.20 (sec) , antiderivative size = 100, normalized size of antiderivative = 1.30 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=-\frac {32 \, {\left (\cos \left (f x + e\right )^{3} - 3 \, \cos \left (f x + e\right )\right )} a^{2} c + 3 \, {\left (12 \, f x + 12 \, e + \sin \left (4 \, f x + 4 \, e\right ) - 8 \, \sin \left (2 \, f x + 2 \, e\right )\right )} a^{2} c - 24 \, {\left (2 \, f x + 2 \, e - \sin \left (2 \, f x + 2 \, e\right )\right )} a^{2} c + 96 \, a^{2} c \cos \left (f x + e\right )}{96 \, f} \]

[In]

integrate(sin(f*x+e)*(a+a*sin(f*x+e))^2*(c-c*sin(f*x+e)),x, algorithm="maxima")

[Out]

-1/96*(32*(cos(f*x + e)^3 - 3*cos(f*x + e))*a^2*c + 3*(12*f*x + 12*e + sin(4*f*x + 4*e) - 8*sin(2*f*x + 2*e))*
a^2*c - 24*(2*f*x + 2*e - sin(2*f*x + 2*e))*a^2*c + 96*a^2*c*cos(f*x + e))/f

Giac [A] (verification not implemented)

none

Time = 0.37 (sec) , antiderivative size = 59, normalized size of antiderivative = 0.77 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=\frac {1}{8} \, a^{2} c x - \frac {a^{2} c \cos \left (3 \, f x + 3 \, e\right )}{12 \, f} - \frac {a^{2} c \cos \left (f x + e\right )}{4 \, f} - \frac {a^{2} c \sin \left (4 \, f x + 4 \, e\right )}{32 \, f} \]

[In]

integrate(sin(f*x+e)*(a+a*sin(f*x+e))^2*(c-c*sin(f*x+e)),x, algorithm="giac")

[Out]

1/8*a^2*c*x - 1/12*a^2*c*cos(3*f*x + 3*e)/f - 1/4*a^2*c*cos(f*x + e)/f - 1/32*a^2*c*sin(4*f*x + 4*e)/f

Mupad [B] (verification not implemented)

Time = 14.06 (sec) , antiderivative size = 250, normalized size of antiderivative = 3.25 \[ \int \sin (e+f x) (a+a \sin (e+f x))^2 (c-c \sin (e+f x)) \, dx=\frac {a^2\,c\,x}{8}-\frac {{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2\,\left (\frac {a^2\,c\,\left (3\,e+3\,f\,x\right )}{6}-\frac {a^2\,c\,\left (12\,e+12\,f\,x-16\right )}{24}\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6\,\left (\frac {a^2\,c\,\left (3\,e+3\,f\,x\right )}{6}-\frac {a^2\,c\,\left (12\,e+12\,f\,x-48\right )}{24}\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4\,\left (\frac {a^2\,c\,\left (3\,e+3\,f\,x\right )}{4}-\frac {a^2\,c\,\left (18\,e+18\,f\,x-48\right )}{24}\right )+\frac {a^2\,c\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}{4}-\frac {7\,a^2\,c\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3}{4}+\frac {7\,a^2\,c\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5}{4}-\frac {a^2\,c\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^7}{4}+\frac {a^2\,c\,\left (3\,e+3\,f\,x\right )}{24}-\frac {a^2\,c\,\left (3\,e+3\,f\,x-16\right )}{24}}{f\,{\left ({\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2+1\right )}^4} \]

[In]

int(sin(e + f*x)*(a + a*sin(e + f*x))^2*(c - c*sin(e + f*x)),x)

[Out]

(a^2*c*x)/8 - (tan(e/2 + (f*x)/2)^2*((a^2*c*(3*e + 3*f*x))/6 - (a^2*c*(12*e + 12*f*x - 16))/24) + tan(e/2 + (f
*x)/2)^6*((a^2*c*(3*e + 3*f*x))/6 - (a^2*c*(12*e + 12*f*x - 48))/24) + tan(e/2 + (f*x)/2)^4*((a^2*c*(3*e + 3*f
*x))/4 - (a^2*c*(18*e + 18*f*x - 48))/24) + (a^2*c*tan(e/2 + (f*x)/2))/4 - (7*a^2*c*tan(e/2 + (f*x)/2)^3)/4 +
(7*a^2*c*tan(e/2 + (f*x)/2)^5)/4 - (a^2*c*tan(e/2 + (f*x)/2)^7)/4 + (a^2*c*(3*e + 3*f*x))/24 - (a^2*c*(3*e + 3
*f*x - 16))/24)/(f*(tan(e/2 + (f*x)/2)^2 + 1)^4)